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Abstract 

We study the recently introduced notion of output-input stability, which is a robust variant of the 



u 
o 

c"| minimum-phase property for general smooth nonlinear control systems. The subject of this paper 

is developing the theory of output-input stability in the multi-input, multi-output setting. We show 
that output-input stability can be viewed as a combination of two system properties, one related to 
detectability and the other to left-invertibility. For systems affine in controls, we provide a necessary 
and sufficient condition for output-input stability, which relies on Hirschorn's nonlinear structure 
algorithm. 

> 1 

in 
o 

(N ■ 1 Introduction 
in 
O 

For systems with inputs, two properties of interest are asymptotic stability under zero inputs and 
bounded state response to bounded inputs. It is well known that for linear time-invariant systems the 
first property implies the second one, but for nonlinear systems this is not the case. The notion of input- 
to-state stability (ISS) introduced in [11] captures both of the above properties; it requires that bounded 



inputs produce bounded states and inputs converging (or equal) to zero produce states converging to 



zero. 



Dual concepts of detectability result if one considers systems with outputs. For linear systems, one of 
equivalent ways to define detectability is to demand that the state converge to zero along every trajectory 



for which the output is identically zero. The notion of output-to-state stability (OSS) introduced in [12] 
is a robust version of the detectability property for nonlinear systems and a dual of ISS; it requires that 
the state be bounded if the output is bounded and converge to zero if the output converges to zero. 

The present line of work is concerned with the minimum-phase property of systems with both inputs 
and outputs. A linear system is minimum-phase if whenever the output is identically zero, both the 
state and the input must converge to zero; in the frequency domain, this is characterized by stability of 
system zeros. Byrnes and Isidori ^ provided an important and natural extension of the minimum-phase 
property to nonlinear systems (affine in controls). According to their definition, a system is minimum- 
phase if its zero dynamics — the internal dynamics of the system under the action of an input that holds 
the output constantly at zero — are asymptotically stable. 

The above remarks suggest that to complete the picture, one should have a robust version of the 
last property, which should ask the state and the input to be bounded when the output is bounded and 
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to become small when the output is small. Such a concept was proposed in the recent paper || under 
the name of output-input stability. In the spirit of ISS, output-input stability requires the state and the 
input of the system to be bounded by a suitable function of the output and derivatives of the output, 
modulo a decaying term depending on initial conditions. This property is in general stronger than the 
minimum-phase property]]] defined in Q . Output-input stability can be studied with the help of the tools 
that have been developed over the years to study ISS, OSS, and related notions (such as Lyapunov- 
like dissipation inequalities); the minimum-phase property, on the other hand, is investigated using 
different techniques (such as computation of zero dynamics and normal forms). This makes output- 
input stability an appropriate alternative notion to use in those situations where the minimum-phase 
property is insufficient or difficult to check. 

The results of || provide a fairly complete theory of output-input stable single-input, single-output 
(SISO) nonlinear control systems. In this paper we continue to study the output-input stability property 
for multi-input, multi-output (MIMO) systems. Our goal is to investigate a connection between output- 
input stability and structural properties of control systems which have been studied in the context of 
system inversion. In particular, we show the relevance of Hirschorn's nonlinear structure algorithm Q 
(which is an extension of Silverman's linear structure algorithm (8|, ||]) in establishing output-input 
stability. After providing necessary definitions and preliminary results in Section |2| and reviewing the 
nonlinear structure algorithm in Section |3|, we derive our main result for affine systems in Section ||| and 
then discuss some extensions in Section |5|. 



2 Background and preliminary results 

Consider the system 

x = f(x, u) 

. (1) 
V = h(x) 

where the state x takes values in M n , the input u takes values in M m , the output y takes values in 
MP (for some positive integers n, m, and p), and the functions / and h are smooth. In this paper we 
restrict admissible input (or "control") signals to be at least continuous. For every initial condition 
x(0) and every input u(-), there is a solution x(-) of ([[]) defined on a maximal interval [0, T max ), and the 
corresponding output y(-). We write C k for the space of k times continuously differentiable functions 
u : [0, oo) — ► M m , where k is some nonnegative integer. Whenever the input u is in C k , the derivatives 
y, y, . . . , exist and are continuous; they are given by 

yW( t )=H i (x(t),u(t),... M'- 1 *®), i = l,...,fc + l, t€[0,T max ) (2) 

where for i = 0, 1, . . . the functions Hi : M n x (M m ) 4 — > MP are defined recursively via Hq := h and 

H i+1 (x,u ,... ,Ui) := ^-f(x,u ) +^2^: u j+i 

X j=o Uj 



1 Strictly speaking, this statement only makes sense for systems affine in controls, because otherwise the minimum-phase 
property is not defined. For example, the scalar system y = 1 + y 2 + u 2 is output-input stable (because |w| < ^/y) but not 
minimum-phase (in fact, no input can hold the output at zero). 
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(here the arguments of Hi are x € M n and uq, . . . , € M m ). Given a nonnegative integer k, we will 
denote by y k the M p ( fc+1 )-valued signal 

/ v \ 

y 

provided that the indicated derivatives exist. We will let || • ||r w denote the supremum norm of a signal 
restricted to an interval [a, b], i.e., ||^||r a h := sup{|z(s)| : a < s < 6}, where | • | is the standard Euclidean 
norm. 

According to Definition 1 of the system (|l|) is called output-input stable if there exist a positive 
integer N, a class ICC function^ /3, and a class /Cqo function 7 such that for every initial state x(0) and 
every input u G C N ~ l the inequality 



x{t) 
u(t) 



</3(|*(o)|,t) +7(l|y iY ll [ o lt] ) 



(3) 



holds for all t in the domain of the corresponding solution. (The assumption that u belongs to C N ^ 
is made to guarantee that is well defined, and can be weakened if the function Hjq is independent 
of un-i-) As discussed in [||, the concept of output-input stability finds applications in feedback 
stabilization, adaptive control, and other areas. 

It is perhaps best to interpret output-input stability as a combination of two separate properties of 
the system. The first one is expressed by the inequality 



\x(t)\ < P(\x(Q)\,t) +7 (||y^|| [0it] ) 



(4) 



and corresponds to detectability (OSS) with respect to the output and its derivatives, uniform over 
inputs; this property was called "weak uniform O-detectability of order N" in ||]. The results of j|, |l2[ 
imply that the system (|l]) is weakly uniformly O-detectable of order N if there exists a continuously 
differentiable, positive definite, radially unbounded function V : W n — > IR that satisfies 



dV 

-r^f{x,uo) < -a(\x\) + x 



H Q (x) 
^H N (x,u , ■ ■ ■ ,u N -i) 



Vx,u , . . . ,u N -! 



for some functions a,x G ^Coo- As discussed in ||, the class of weakly uniformly O-detectable systems 
includes all affine systems in global normal form with ISS internal dynamics. 

The second ingredient of the output-input stability property is described by the inequality 



\u(t)\<[3(\x(0)\,t)+ 7 (\\y N \\ m ) 

which says that the input should become small if the output and its derivatives are small, 
speaking, this suggests that the system has a stable left inverse in the input-output sense, 
uniform detectability, this property does not seem to admit a Lyapunov-like characterization. 
SISO case it is closely related to the existence of a relative degree; see [||, Theorem 1]. In 
however, this second property needs to be understood better, which is precisely the goal of the 



2 Recall that a function a : [0,oo) — > [0, 00) is said to be of class K. if it is continuous, strictly increasing, and 
If a £ K is unbounded, then it is said to be of class /Coo. A function f3 : [0, 00) x [0, 00) — + [0, 00) is said to be of 
if f3(-,t) is of class K, for each fixed t > and /3(s, t) decreases to as t — *■ 00 for each fixed s > 0. 



(5) 

Loosely 
Unlike 
In the 
general, 
present 

a(0) = 0. 
class K.C 
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paper. In what follows, we formulate and study a useful property which, in combination with (||), 
yields (|). 

Let us say that the system (|l]) has the input- bounding property if there exist a positive integer k* 
and two class /Cqo functions p\ and pi such that we have 



\uo\ < Pi(\x\) + Pi 




Vx,uo,...,Uk*-i. (6) 



(We could omit the first component Hq(x) of the argument of p 2 without loss of generality, but will 
keep it for notational convenience.) Defined in this way, the input-bounding property represents a 
functional relation between the input and state variables. The next result recasts this property in terms 
of trajectories of the system. 

Lemma 1 The system (|l|) has the input-bounding property if and only if there exist a positive integer 
k* and two class /Coo functions p\ and p% such that for every initial condition and every input u G C fc * _1 
the inequality 

\u(t)\ < Pl (\x(t)\) + p 2 (\y k * (t)\) (7) 
holds for all t in the domain of the corresponding solution. 

Proof. In view of (|2|), it is clear that (|6|) implies (0), with the same k*, pi, and p%- To show the 
converse, suppose that (0) is violated for some x,uq, . . . , itfc*_i. Take x to be the initial condition and 
apply an input u satisfying ?/W(0) = m, i = 0, . . . , k* — 1. Then it is easy to see that (0) does not hold 
for small t. □ 

We point out that the input-bounding property resembles in its appearance the notion of relative 
degree as defined in ||], but is actually much less restrictive (especially for MIMO systems). The 
next result reveals the connection between output-input stability, weak uniform O-detectability, and the 
input-bounding property; it was essentially proved in @ . 

Lemma 2 Suppose that the system (Q) is weakly uniformly 0-detectable of order k: 

|x(t)|<^|x(0)|,t)+7(||y'|| M ) (8) 

where k > 0, (3 G ICC, and 7 G /Cqo. Suppose also that it has the input-bounding property as defined 
above. Then it is output-input stable, with N = m&x{k, k*}. 



Proof. Since (|l|) has the input-bounding property, by Lemma [T| the inequality (|7|) holds with pi,p2 G 
/Cqq. Combining this with (||) and using the simple fact that for every class /C function p and arbitrary 
numbers si,S2 > one has p(s\ + S2) < p(2si) + ,0(252), we arrive at the inequality (||) with N := 
max{&, &*}, (3(s,t) := pi(2/3(s, t))+/3(s, t), and7(s) := pi(2 ; y(s))+p2(s)+ ; y(s). Thus (jl]) is output-input 
stable. □ 

The above observation explains the importance of the input-bounding property. As we will see 
later, a natural way of checking this property for systems affine in controls is provided by Hirschorn's 
nonlinear structure algorithm. 
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3 Nonlinear structure algorithm 



In this section and the next one we assume that m < p and the system (|l|) takes the form 

x = fix) + G(x)u 

y = h{x) 

Its dynamics can also be written in more detail as 

rn 

x = fix) + ^2gi(x)ui. 

i=l 

We assume that /(0) = and h(Q) = (although the second assumption is only made for convenience 
and can be removed). For / a vector field and R a scalar- valued function, LfR denotes the derivative 
of R along /. We will also let Lf act on vectors and matrices componentwise; for example, if R is 
matrix- valued, then LfR is a matrix- valued function whose ijth component is given by LARij). All 
functions are assumed to have the smoothness required for all relevant derivatives to exist. Dimensions 
of vectors and matrices will be omitted when clear from the context. 

To make the paper self-contained, we now describe (a global version of) Hirschorn's structure algo- 
rithm from Q. It iteratively constructs a sequence of vectors z k £ W which take the form 

z k = M k (x)y k 

(i.e., their components are x-dependent linear combinations of the original outputs and their derivatives) 
and satisfy 

z k = h k {x) + J k {x)u 

where 



J k (x) 







with J k (x) an r k x m matrix of constant rank r k > 0. Partitioning all vectors accordingly, we can 
rewrite this as 



z k \ = (Mk(x)\ , = (h{x)\ (3 k ix 



V \M k (x)) y = UIwJ + { ) < 10 > 

The initial data for k = are given by Mq = I, zo = 0, zq = y, Jiq = 0, ho = h, Jq = 0, and ro = 0. 
To obtain z k+ \, we differentiate z k along solutions of the system. On one hand, we have 

d m 
k k = —(M k (x)y k ) = L f M k (x)y k + ^ L gi M k {x) Ui y k + M k {x)y k . 

i=l 

Assumption 1. L g M k (x) = for all i 

If this assumption holds, then the dependence on u disappears and we obtain 

k k = L f M k (x)y k + M k (x)y k . (11) 
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On the other hand, 



dt 

where 



d 

Zk = -jj.h k {x) = Lfh k (x) + S ^Lg i h k (x)u i = L f h k (x) + J k (x)u 

i=l 



J k (x) := ^G(x). 



~z k \ _ ( h k (x) \ + / J k (x) ^ (12) 



We thus arrive at the equation 

\L f h k (x)J ' \J k (x) 
/J k (x)\ 

Assumption 2. The matrix I ~ I has a constant rank r k +i and a fixed set of r k+ \ rows that are 



Jk(x] 

linearly independent for all x. 

This assumption implies that there exists a p x p permutation matrix E k (independent of x) such 
that 

'Jk(x)\ = /J fc+1 (x) N 
Jk(x)J \J k +i(x) 



E k 



where J k+ i(x) is an r k+ \ x m matrix of rank r k+ \. (In fact, typically it suffices to reorder the rows of 
J k (x), since the rows of J k (x) are already linearly independent.) Therefore, a (p — r k +i) x r k+ \ matrix 
F k (x) can be found such that F k (x)J k +i(x) + J k +\(x) = 0. Consider the p x p matrix 

Multiplying both sides of (Il2|) by R k (x), we obtain 



^)if;)-^)(4t(i)) + ( Jt+ o w )«- 



The iteration step is completed by setting 



z k+ i := R k {x) (l k ) , h k+ i(x) := R k (x) ( J 1 ^) , 



and noticing that from ( |lO|) and ( |TT|) we have = M k+ \ (x)y k+l where 

/ M fc (x) \ 

Mk ^ x) =\( Lf M k {x) o) + (o M fc (x; 

It is easy to see that r k+ \ > r k for all k. If for some k we have r k = m, then the algorithm terminates. 
We denote the smallest such k by k* or let k* = oo if the above condition is never met. 

It is not hard to show that Assumption 1 is fulfilled if L gi L^F k {x) = for 0<j + k<k* — 2 and 
1 < i < m Theorem 2]. (In particular, it is trivially satisfied if k* = 1.) Hirschorn also gives sufficient 
conditions for this in terms of the original system [§, Corollary 1]. Assumption 2 represents a global 
regularity condition. We discuss possible ways of relaxing these assumptions in Section The following 
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example illustrates the application of the above algorithm to a system satisfying both assumptions at 
each step. 



Example 1. Consider the system 



X\ = U\ 

X 2 = ^3 + x^u\ 



X3 = u 2 (13) 
X4. = — X4. H- x-^ 
y = (xi,x 2 ) T 



We have 



yi ) = r)+ 1 °) r 

m \x?j \x 4 01 \u 2 



yi \ A (A /«! 

x 4 - x? - x 4 y\ + y 2 ) \0 l) \u 2 



so ri = 1. Multiplying both sides by Rq = ( ^ ) yields 

y— X4 ly 

-xm + mj \X3/ V° °/ W 

Differentiating the bottom component, we obtain 

(14) 

hence r 2 = 2 and the algorithm terminates with k* = 2. □ 

4 Main result 

We are now ready to state and prove our main result for affine systems. 

Theorem 1 Let Assumptions 1 and 2 hold for each k > 0. Then the system (^) is output-input stable 
if and only if it is weakly uniformly 0-detectable (of some order k) and the algorithm of Section^ gives 
k* < 00. 

Proof. We first prove sufficiency. If k* < 00, then the system is left-invertible, i.e., one can explicitly 
solve for u as a function of x and y fc * . Indeed, we have 

Mk*{x)y k * = h k *{x) + J k *(x)u 

where J k *{x) is square and invertible. Multiplying both sides by J^}{x), we obtain 

u = -J^}(x)h k *(x) + J^}{x)M k ,(x) Y k ' =: A(x) + B(x)y k * . 

It is easy to show that /(0) = implies ^4(0) = 0. It follows that 

L4(x)|< 7l (M), \B(x)\<\B(0)\+<n(\x\) 
for some 71 , 72 £ /Coo • This yields 

M < 7 i(M) + (|B(o)[ + 72(H)) |y fc *| < n(\x\) + \B(o)\\y h *\ + ^(M)) 2 + ^|y fc T 
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and so the inequality ([?]) holds with pi(r) := 71 (r) + \{^2{f)) 2 and p 2 {r) := |5(0)|r + \r 2 . By Lemma [l] 
the system has the input-bounding property, and Lemma |2| guarantees output-input stability with 
N = max{A;, k*}. 

Necessity follows from the correspondence between Hirschorn's structure algorithm and the zero 
dynamics algorithm, which is explained, e.g., in |?], Sections 11.1-11.2] or ||, Section 6.1]. If k* = 00, 
then in some neighborhood of the origin there exist a vector a(x) € M m and an m x I matrix (3(x), 
where < I < m, such that every feedback law of the form u = a(x) + (3(x)v, v € M! holds the output 
constantly at zero (for a proper choice of initial conditions). Clearly, this violates the inequality (||), 
hence the system cannot be output-input stable. The fact that output-input stability implies weak 
uniform O-detectability is an immediate consequence of the definitions. □ 

Example 1 (continued). Consider again the system (|l3|). From the formula ( |l~4|) we conclude that 
< \yi\ and 

\u 2 \ < -{xa - x x ) + -y x + -x 4 + -y x + \y 2 \ 

hence the system has the input-bounding property. It is also weakly uniformly O-detectable of order 1, 
as is seen from the bound |x3| = \ij2 — £43/1 1 < \y% \ + |a;4||yi| and the fact that the equation for x 4 is ISS 
with respect to x\. Therefore, the system ( |l3D is output-input stable. □ 

A sufficient condition for termination of Hirschorn's algorithm, and consequently for the input- 
bounding property, is the existence of a uniform (vector) relative degree as defined, e.g., in Q. This 
condition is by no means necessary; note that the system considered in the above example does not 
have a uniform relative degree. 

The next example illustrates what can happen when Assumption 2 is violated. 

Example 2. Consider the system 

X\ = U\ 

x 2 =xz + x 2 u 2 

X3 = u 2 (15) 

y = (x 1 ,x 2 ) T 



We have 



* i ) = (°) + ( 1 r 

2/2/ \X3/ V° X 2/ \U2 



and the rank of the matrix on the right-hand side drops from 2 to 1 when x 2 = 0. It is not difficult to 
show that we can pick a bounded sequence of initial states along which x 2 (0) converges to 0, a sequence 
of values of 1x2(0) converging to 00, and appropriately chosen sequences of values for u 2 (0), ^2(0), ... 
such that the derivatives y 2 (0), 2/2(0), ... are all kept at zero. Also, let u\ = so that y\ = 0. This 
implies that the inequality (|5|) is violated for small t, hence the system ( |l~5|) is not output-input stable. 
(The proof of Theorem 1 in || contains a general argument along these lines.) □ 

It is instructive to note that both the system ( |l~3|) considered in Example 1 and the system (|l5|) 
considered in Example 2 are minimum-phase, with zero dynamics in both cases being given by xa = —xa- 
An important fact not elucidated by the zero dynamics is that the minimum-phase property of the 



system (13) is "robust" (small y, y, ... force x and u to be small) while the minimum-phase property 



of the system ( |15|) is "fragile" (small y, y, ... can correspond to arbitrarily large u). 
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5 Some extensions 



Assumptions 1 and 2 are automatically satisfied for SISO affine systems with uniform relative degree 
in the sense of ||. In the SISO case, output-input stability actually implies the existence of relative 
degree for a class of systems which includes systems affine in controls; see Theorem 1]. We will see 
below that in contrast with that result, for MIMO systems Assumptions 1 and 2 are not necessary for 
output-input stability. 

Suppose that at some step k Assumption 1 does not hold. Then we can use the following modification 



of Hirschorn's algorithm, proposed in [1C]. Let us collect all u-dependent terms on the right-hand side: 



L f M k {x)y k + M k {x)y k = L f h % (x) + ^L ft %(x)u; - ^ L g .M- k (x)u iy k =: L f h k (x) + J k (x,y k )u. 



»=i 



Then we can follow the same steps as before and define matrices R/~, k > k and k > k which will 
depend rationally on y k , and vectors z k , k > k which will depend rationally but not necessarily linearly 
on y k . Assuming that Assumption 2 is still satisfied with "for all x,y fc " replacing "for all x" , and that 
k* < oo, we can bound u in terms of x and a rational function of y k * (cf. M, Example 11.20]). This in 
general does not give the input-bounding property. However, in many cases the above dependance on 
y fc * is actually polynomial, and then the input-bounding property does hold. The next example serves 
to support this observation. 

Example 3. Consider the system 

±i = U\ 
X 2 = ^3 + X 2 U\ 
%3 = U 2 



(16) 



y = (x 1 ,x 2 ) 



We have 



+ 



so r\ = 1. Multiplying both sides by Rq 



yi 

-%2yi + m 



o 

1 

-X 2 1 



X3 



1 

X 2 



gives 



+ 



1 




Differentiating the bottom component, we obtain 



yi 



1 o 
o 1 



-x 3 - x 2 ui)yi - x 2 y\ + y 2/ 
Moving the tii-dependent term from the left-hand side to the right, we can rewrite this as 



2/1 

-x 3 yi - x 2 y\ + y 2 
-l 



Finally, multiplying both sides by 



1 







x 2 yi 1 



1 

x 2 yi 1 





1 



-x 2 yi 1 



y\ 



-x 2 yf - x 3 yi - x 2 y\ + y 2 



, we arrive at 
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and now the input-bounding property can be shown as in Example 1. Weak uniform O-detectability can 
also be verified along the same lines as in Example 1. Therefore, the system (|jl|) is output-input stable. 
We leave it to the reader to check that if we add another state to the system by replacing the equation 
x 3 = U2 with ±3 = X4 and adding the equation £4 = U2, then the conclusion remains unchanged. □ 

Note that Assumption 1 was not used in the proof of the necessity part of Theorem |]. It is not hard 
to see that necessity is preserved if this assumption is removed and the structure algorithm is modified 
as described above. 

We now demonstrate that a system can be output-input stable even if Assumption 2 is not satisfied. 
(This is to be contrasted with Example 2.) 

Example 4. The system 



Xl 


= Ul 


X 2 


= X 5 + X4U2 


X 3 


= X4 




= U 2 


x 5 


= u 3 


y 


= (xi,x 2 ,x 3 ) T 




is output-input stable, as can be seen from the formulas u% = yi, u% = j/3, U3 = j/2 — y| — 2/3 2/3, X4 = y 3 , 
and X5 = 2/2 — 2/32/3- However, when we try to apply Hirschorn's algorithm, we obtain 

fi o\ /V 

+ x 4 1*2 
\0 0/ \u 3/ 

and the matrix on the right-hand side does not have a constant rank0. □ 

It is straightforward to obtain a local variant of Theorem |]. Namely, if Assumption 2 holds for all k 
only in a neighborhood of a given equilibrium — say, the origin — and if in the definition of output-input 
stability we only require the inequality (|J) to be satisfied for initial conditions close to the origin and for 
small t, then the theorem is still valid. (Note also that locally, the second statement in Assumption 2 
automatically follows from the first.) For analytic systems, Assumption 2 is guaranteed to hold on an 
open and dense subset of M. n (the complement of a singular set). Local versions of output-input stability 
may be useful for control design purposes away from singularities (cf. 

The above results can be used to establish output-input stability of some nonaffme systems. Note 
that to have the input-bounding property, we only need to be able to bound — and not necessarily solve 
for — the input in terms of the state and derivatives of the output. As a simple generalization, consider 
a system of the form 



x 



f( x ) + ^Z9i{x)li{ui) 



where the functions ji, i = 1, . . . , m are bounded from below by some class /Coo functions. It is easy to 
show that if the associated "virtual input" system 



m 

x = f{x) + ^gi{x)vi 
1=1 



This system has trivial zero dynamics (x = 0) but x = is not a regular point of the zero dynamics algorithm. 
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is covered by the sufficiency part of Theorem |l], then the original nonaffine system is output-input stable. 
Of course, left-invertibility of the virtual input system is not necessary. One can even have more inputs 
than outputs; for example, the scalar system y = u\ + u\ is clearly output-input stable. 

There has been a lot of activity in the 1980s on establishing invertibility conditions for nonlinear 
systems in terms of differential geometric concepts — particularly, controlled invariant submanifolds. 
These and related developments are documented in |3|, Chapter 6], J/], Chapter 11], and the references 
therein, and build on earlier results from geometric theory of linear systems ||. This body of work 
helps one to better understand the above constructions and reveals their coordinate independence (see, 
e.g., the proof of the necessity part of Theorem |l|). However, these results seem to be of limited utility 
in establishing output-input stability, as they do not capture properties such as the one expressed by 
Assumption 1 or the polynomial dependence discussed earlier in this section. It would be interesting to 
gain more insight into the questions studied here with the help of geometric tools. 
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